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Abstract 



We present a measurement of the like-sign dimuon asymmetry in semilep- 
tonic 5-hadron decays, performed using 6.1 fb _1 of pp collisions recorded 
with the DO detector at a center-of-mass energy y/s = 1.96 TeV at the 
Fermilab Tevatron collider. This measured value is A sl = [—0.957 ± 
0.251 (stat) ± 0.146 (syst)] %, which disagrees with the Standard Model 
prediction at a statistical level of 3.2 a, and provides the first evidence of 
anomalous CP violation in the mixing of neutral B mesons. 
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1 Introduction 



During the big-bang origin of the universe, fundamental symmetry principles imply 
that particles and antiparticles were created in equal quantities. However, current 
observations show the universe to be matter dominated. Baryogenesis provides a 
mechanism to generate the observed asymmetry and requires a number of conditions, 
one of which is the violation of CP symmetry pp. The Standard Model (SM) naturally 
includes CP violation (CPV) through complex phases induced by quark mixing. 
While measurements in the K° and B d systems show excellent agreement with the 
SM predictions [2], the observed level of CP violation is far too small to yield the 
observed matter-antimatter asymmetry [3]. As such it is important to search for 
new sources of CPV beyond the SM. The B° s system is one possible source, since 
the effects of CPV are expected to be small in the SM [4], and can be significantly 
enhanced by new physics models [5j El [TJ El |9] . 

CPV effects can be induced via the mixing process whereby a neutral B meson 
oscillates into its antiparticle B. In this case, the experimental observable is the 
asymmetry between the processes B° s d — > B s d — >■ / and B sd ^B® d ^-f, where / is 
the final state decay product. In particular, for semileptonic decays to a muonic final 
state, the flavor specific asymmetry is defined as: 

q _ T (B q -^B q ^ fx+X) - T(B q ^B q ^ fj,-X) 



sl T(B q ^B q ^ fi+X) + T(B q ^B q ^ fi-X) ' 

For the case where the original flavor (B d , B®) is not distinguished, the asymmetry 
has contributions from both systems, and is denoted a b sl = 0d a ti + Ps^i- At the 
Tevatron accelerator, the coefficients have been measured as (3d = 0.506 ± 0.043 and 
P. s = 0.494 ± 0.043 [2]. The asymmetries are related to the meson mixing parameters 
by the relation: 

Ar 

< = XW q ^^ (!) 

where <p q is the CP-violating phase, and AM q (AT g ) is the mass (width) difference 
between the eigenstates of the neutral B q meson systems. The SM predicts af { = 
(-4.8+};°) x 10" 4 , a s s ! = (2.06 ± 0.57) x 10~ 5 [4J, and hence: 

4(SM) = [-0.0231°;^] %• (2) 

Experimentally, it is a challenge to measure such an asymmetry to high precision, 
in part because of dilution from background processes where the mesons don't mix 
prior to decay. A more sensitive parameter is the like-sign dimuon asymmetry from 
semileptonic decays: 

A b = Nb++ ~ ^ (3) 
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where N b are the number of events containing two 6-hadrons decaying semilep- 

tonically into two positive (negative) muons. For an original bb state to produce 
like-sign dimuons in the final state, one of the quarks must have hadronised to a 
neutral meson, and then mixed prior to decay, while the other must decay without 
mixing. In fact, there is a second source of such muons: sequential semileptonic 
decays b —> c\i~X — > fi + fi + X'. Such events are suppressed in this analysis by the 
requirement M(/z/i) > 2.8 GeV, and their contributions are taken into account. The 
fundamental asymmetry being probed is unchanged by requiring a second muon, i.e. 

^si — a si- 

In this document, we present a measurement of A h sll which takes advantage of the 
two different ways it can be extracted to minimise the total experimental uncertainty 
on the asymmetry. This study was performed using 6.1 fb _1 of data recorded with the 
DO detector [10J at the Fermilab Tevatron proton-antiproton (pp) collider, operating 
at a center-of-mass energy of 1.96 TeV. More detailed information can be found in 
the official publication jllj . 

The DO experiment is well suited to the investigation of the small effects of CP 
violation because the periodic reversal of the DO solenoid and toroid magnetic field 
polarities results in a cancellation of most detector-related asymmetries. In addition, 
the pp initial state is a CP eigenstate, and the high center-of-mass energy provides 
access to mass states beyond the reach of the 5-factories running at y/s = M(T(4S')). 



2 Measurement Overview 

To extract A^ h the raw asymmetries must first be measured. The like-sign dimuon 
charge asymmetry A is defined as 

A = N++ + N- > (4) 

where N ++ ( > is the number of events containing two positive (negative) muons 
which fulfill the kinematic requirements (see below). This quantity is related to 
by: 

A = K-Al l + A hkg . (5) 

Here A^ g represents the contribution to the asymmetry from background processes, 
defined as those in which the muon arises from the decay of a long-lived particle, and 
where detector asymmetries can therefore play a role. The coefficient K accounts for 
dilution from signal-like events (i.e. weak decays of heavy quarks) without B® mixing, 
which only contribute to the denominator in Eq. pi) . 
Similarly, the inclusive muon asymmetry is defined as 

n+-n- , , 



n + + n~ 



k ■ A sl + a bkg , (6) 
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with the lower-case symbols representing the equivalent inclusive muon parameters 
as those described for the dimuon case above. 

The general strategy of the measurement is as follows. The raw asymmetries A and 
a are first obtained by event counting, following application of all selection criteria. 
The parameters A^g, a^kg, K and k are determined in bins of muon transverse 
momenta, using data-driven methods wherever possible. Equations (TJH6]) are then 
used to extract A^ separately for the inclusive and dimuon cases. Finally, the two 
measurements are combined, taking advantage of correlations in their uncertainties, 
to yield a single measurement of with significantly improved precision. 

We define all muons from weak decays of b and c quarks as signal, and all other 
muons (dominated by pion and kaon decay-in-flight) as background. The dilution co- 
efficients are determined from Monte Carlo simulation, using the branching fractions 
and momentum spectra of signal processes to evaluate their relative contributions 
to the signal sample. Apart from the asymmetry-inducing semileptonic B decays, 
contributions from the following processes are taken into account in the simulation: 
sequential decays b — > c — > fi + X, decays b — > ccq, cc and bbcc production, and the 
dimuon decays of narrow resonances, such as rj, u, p°, 0(1020) and J/"0. The obtained 
values are k = 0.070 ± 0.006 and K = 0.486 ± 0.032, showing that the raw dimuon 
asymmetry is more sensitive to A^ than the inclusive muon asymmetry, as expected 
from the reduction in background contributions. 

The largest background asymmetry arises from the fact that K + and K~ mesons 
interact differently with the material of the detector (which is clearly matter-antimatter 
asymmetric), and thus their decay rates into positive and negative muons are not 
identical. Background contributions from kaon and pion decay- in-flight, hadronic 
punchthrough to the muon system, and residual reconstruction asymmetries, are all 
measured using collision data. 

3 Event Selection 

The inclusive muon and like-sign dimuon samples are collected from events satisfying 
single and dimuon triggers, respectively. Muon candidates are defined by matching 
central tracks with a 'local muon' track reconstructed in the muon system and by 
applying tight quality requirements aimed at reducing backgrounds such as cosmic 
rays, beam halo, and mis-matched muon-track combinations. All muons must have 
transverse momentum satisfying 1.5 < pr < 25 GeV, and lie in the geometrical 
region corresponding to pseudorapidity \r]\ < 2.2. For muons with px < 4.2 GeV, 
an additional requirement is applied to the longitudinal component \p z \ > 6.4 GeV, 
to ensure that the muon has sufficient penetrating power to reach the outermost 
layers of the muon system. The transverse impact parameter of the muon track 
relative to the reconstructed primary pp interaction vertex (PV) must be smaller than 
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0.3 cm, and this point-of-closest-approach must lie within 0.5 cm of the PV along 
the beam axis. Muons fulfilling the above requirements define the inclusive muon 
sample, used to determine a in Eq. [6j with no constraint on the number of muons per 
event. The dimuon sample contains all events with at least two muon candidates of 
the same charge. To reduce the contribution of background from sequential decays 
B — > jiDX — > /i/iX', the two muons must have an invariant mass above 2.8 GeV. 

4 Detector Asymmetries 

The background asymmetries a^ g and A^ g comprise contributions from several 
sources. These include secondary muons from decays-in-flight of charged kaons and 
pions; fake muon signals from hadronic 'punch-through' of kaons, pions, and protons 
to the muon system; and residual detector asymmetries. For example, the inclusive 
muon background asymmetry is expressed as: 



Here, the first three terms account for asymmetries from kaons, pions, and protons, 
respectively. In each case, the total asymmetry is the product of the fraction of 
muons arising from this source, and the asymmetry a« associated with that source. 
The last term accounts for any remaining detector asymmetry, 5, on the signal muon 
sample. The dominant source of asymmetry is the kaon term, since the K~ meson 
has a significantly shorter interaction length in matter than the K + meson [2], and 
therefore has less chance of decaying to a muon. The resulting contribution to the 
raw asymmetry a is therefore positive, and at least a factor ten larger than any of 
the other three terms in Eq. ([7]). 

The asymmetry from kaons is measured in data using K*° — > K^n^ decays. 
The K*° invariant mass distribution is constructed by assigning the charged kaon 
mass to all muons in the inclusive muon sample and combining with another track 
(assigned the pion mass). The difference of distributions for events with K + — > fi + 
and K~ — > yT tracks is shown in Fig. [U demonstrating the large excess of positive 
K — > fi candidates. Simulations are used to correct for the fraction of kaons (~ 6%) 
that decay prior to being reconstructed [11]. The measurement is repeated in an 
independent sample of kaons from <fi — > K + K~ decays, and the results are fully 
consistent. Similarly, the asymmetry for pion or proton tracks misidentified as muons 
(a n , a p ) are measured using samples of K s — > 7r + 7r~ and A — > pn~ decays, respectively. 

The fractions of muons from kaons, fx, is extracted using a similar method, under 
certain theoretical assumptions which are validated in data. Additional information 
from both data and simulation is then used to relate the kaon fraction to the corre- 
sponding pion and proton fractions and f p . 



a bkg — fK ■ a K + fn ■ a TT + fp ■ dp 




(7) 
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After measuring the fractions of muons originating from kaons, pions and protons, 
the remaining signal (heavy flavor) muons constitute [58.1 ±1.4 (stat)±3.9 (syst)]% of 
the inclusive muon sample. For the like-sign dimuon case, a fraction [66.5±1.6 (stat)± 
3.3 (syst)]% of the candiate events have both muons produced from heavy flavor 
decays. 

The effects of detector asymmetries on this measurement are significantly sup- 
pressed by the regular reversal of the solenoid and toroid polarities in the DO detector. 
Using J ftp — > decays, the residual reconstruction asymmetry is measured as a 

function of muon p?, and found to be of order 10~ 3 . More detail on such effects can 
be found in Ref. p2]. 

5 Results 

The uncorrected asymmetries a and A are obtained by counting events in the inclusive 
muon sample (containing 1.495 x 10 9 muons) and like-sign dimuon sample (containing 
3.731 x 10 6 events), respectively, and measured to be: 

a = [+0.955 ± 0.003 (stat)] %, (8) 
A = [+0.564 ± 0.053 (stat)] %. (9) 

After correcting for background (a^, A bkg ) and for the dilutions (k, K) of the signal 
component, we obtain: 

A\ x = [+0.94 ± 1.12 (stat) ± 2.14 (syst)] % (10) 

xlO 2 
% 3000 

§ 2000 
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| 1000 
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M(Kjc) [GeV] 

Figure 1: The difference in the number of events for the K + n~ and K~n + mass 
distributions of K*° candidates in the inclusive muon sample. The solid line represents 
the result of the fit, while the dashed line shows the background contribution. 
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Table 1: Sources of uncertainty on A b y in Eqs. (I10|lll|ll3p . The first eight rows corre- 
spond to statistical uncertainties and the next three rows to systematic uncertainties. 



Source 










m 




m 


A or a (stat) 


0.00066 


0.00159 


0.00179 


K fraction (stat) 


0.00222 


0.00123 


0.00140 


7r fraction 


0.00234 


0.00038 


0.00010 


p fraction 


0.00301 


0.00044 


0.00011 


K asym. 


0.00410 


0.00076 


0.00061 


7r asym. 


0.00699 


0.00086 


0.00035 


p asym. 


0.00478 


0.00054 


0.00001 


detector asym. 


0.00405 


0.00105 


0.00077 


K fraction 


0.02137 


0.00300 


0.00128 


7T, K, p 


0.00098 


0.00025 


0.00018 


multiplicity 


A b Bl dilution 


0.00080 


0.00046 


0.00068 


Total statistical 


0.01118 


0.00266 


0.00251 


Total systematic 


0.02140 


0.00305 


0.00146 


Total 


0.02415 


0.00405 


0.00290 



from the inclusive muon sample and 

A b sl = [-0.736 ± 0.266 (stat) ± 0.305 (syst)] % (11) 

from the like-sign dimuon sample. The uncertainties on the latter measurement are 
much smaller, due to the difference in the dilution coefficients k <C K in Eqs. fl5H6]). 
Since the same background processes contribute to the uncorrected asymmetries a 
and A, their uncertainties are strongly correlated. The systematic uncertainties can 
therefore be significantly reduced by suitable linear combination of the raw asymme- 
tries: 

A' = A — aa. (12) 

The total uncertainty is minimised by scanning the coefficient a over a wide range, 
yielding the optimal value a = 0.959. The final result for the asymmetry A b Y is then: 

A b sl = [-0.957 ± 0.251 (stat) ± 0.146 (syst)] %. (13) 

It differs by 3.2 standard deviations from the SM prediction for A b \ of Eq. fl5]). The 
contributions to the total uncertainty of A^ in Eqs. (I10|lll|ll3p are listed in Table HJ 
showing that the combination reduces the systematic effects to the point where the 
precision is dominated by statistical uncertainties. 
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Figure 2: The observed and expected like-sign dimuon charge asymmetry A as a 
function of dimuon invariant mass. The expected asymmetry is shown for ^4^=0 and 
^=-0.00957. 



6 Cross-Checks 

To test the stability of the measurement under a range of different event selection 
criteria, several consistency checks are performed. In each case, the procedure to 
extract is repeated over a sub-sample of data, defined by applying new criteria 
on detector era, muon and track quality requirements, impact parameter, transverse 
momentum, polar angle and pseudorapidity. While the raw asymmetries a and A vary 
significantly, by up to 140% (as a result of different sample composition), the obtained 
variations of A^ are consistent with purely statistical fluctuations with respect to the 
nominal value. 

The small value of the dilution coefficient k in Eq. (jSJ) shows that the raw inclusive 
muon asymmetry is dominated by background contributions. This is verified by com- 
paring the observed asymmetry with the background expectation a^kg as a function 
of muon pt, and excellent agreement is found. In addition, the expected invariant 
mass dependence of the like-sign dimuon asymmetry is reproduced in data, provided 
that A b sl is fixed to its measured value. The alternative hypothesis A b sl = leads to 
significant discrepancy, as shown in Fig. |2j 
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7 Interpretations and Conclusions 



The asymmetry is a linear combination of the flavor-specific asymmetries a d sl and 
a s sl , and therefore places constraints on the allowed phase space in the a^-a^ plane, 
as shown in Fig. |3j Shown on the same axis are independent constraints on afj (from 
the .B-factories [131 HH US]) and on (from the DO measurement in — > D s [iX 
decays [IS]). Using additional constraints on a d sl allows the CP- violating phase <p s of 
the B° s mixing matrix to be constrained, as shown in Eq.( [T|). This extraction, and 
comparisons with related measurements, are given in Ref. [11] . 

In conclusion, using a sample of pp collison data, corresponding to 6.1 fb _1 of 
integrated luminosity, we have measured the like-sign dimuon charge asymmetry A b sl 
of semileptonic 6-hadron decays: 

A b Bl = [-0.957 ± 0.251 (stat) ± 0.146 (syst)] %. (14) 




-0.04-0.03-0.02-0.01 0.01 

a sl 

Figure 3: Comparison of A^ in data with the SM prediction for a d { and aj. 
Also shown are other measurements of a d Y = —0.0047 ± 0.0046 [T3~| [HI [15] and 
= —0.0017 ± 0.0091 jT6j. The bands represent the ±1 standard deviation un- 
certainties on each measurement. 
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This measurement is consistent with, and supersedes, our previous measurement [12] 
obtained with 1 fb _1 . The asymmetry disagrees with the prediction of the SM by 3.2 
standard deviations, and gives the first evidence for anomalous CP violation in the 
mixing of neutral B mesons. 

8 Future Plans 

The DO collaboration are working on an updated measurement of the asymmetry A^, 
taking advantage of an increased dataset and several improvements to the analysis. 
These are briefly described in this section. 

The largest contribution to the measured uncertainty on A^ is due to the limited 
size of the data sample. The DO detector continues to collect data, with approximately 

9 fb _1 available for analysis at the time of writing. This represents an increase in 
integrated luminosity of ~50% with respect to the published result, therefore an 
updated measurement using this additional data will significantly increase sensitivity 
to non-SM effects, even with no analysis improvements. 

The event selection is being revisited, and a few minor changes have been proposed 
to optimise the signal sample, while maintaining a small background contamination. 
A new method (the so-called 'null fit' descibed in Appendix E of Ref. [IE]) has been 
developed to measure the kaon fractions with greater precision, reducing sensitivity 
to the detector mass resolution and peaking backgrounds in the K*° invariant mass 
distributions. This new method also enables a second channel to be included in the 
Fk/ fk measurement, with consistent results obtained for both. 

Finally, muon impact parameter (IP) information is being used to help identify 
the likely source of the observed asymmetry. For muons with small IP (less than 
100/xm, for instance), the probability that they arise from a B® meson which has 
mixed is small, due to the low mixing frequency associated with the mass difference 
AMd fa 0.5 ps _1 . In contrast, the contribution of muons from decays of B® mesons 
which have mixed prior to decay is relatively stable with IP, due to the much higher 
mixing frequency associated with the mass difference AM S pa 17 ps -1 . As such, by 
dividing the sample into two parts corresponding to small and large impact parameter 
muons, the resulting independent measurements of A^ yield two bands in the a^-a^ 
plane with different slopes. The overlapping region of these two slopes can then be 
directly compared to the SM predictions for both and aj. 

The anticipated timescale for publication is summer 2011. 

References 

[1] A.D. Sakharov, Pisma Zh. Eksp. Teor. Fiz. 5, 32 (1967) [Sov. Phys. JETP Lett. 
5,24 (1967)]. 



9 



[2] C. Amsler et al, Phys. Lett. B 667, 1 (2008). 

[3] M.B. Gavela, P. Hernandez, J. Orloff, O. Pene, Mod. Phys. Lett. A 9, 795 (1994); 
M.B. Gavela, P. Hernandez, J. Orloff, O. Pene, C. Quimbay, Nucl. Phys. B 430, 
382 (1994); P. Huet and E. Sather, Phys. Rev. D 51, 379 (1995). 

[4] A. Lenz and U. Nierste, J. High Energy Phys. 0706, 072 (2007). 

[5] L. Randall and S. Su, Nucl. Phys. B 540, 37 (1999). 



[6] J.L. Hewett, |arXiv:hep-ph/ 9803370 (1998). 
[7] G.W.S. Hou. larXiv:0810.3396l [hep-ph] (2008). 

[8] A. Soni et al, Phys. Lett. B 683, 302 (2010); A. Soni et al, larXiv:1002.0"595l 
(2010) [hep-ph] and references therein. 

[9] M. Blanke, A. J. Buras, A. Poschenrieder, C. Tarantino, S. Uhlig and A. Weiler, 
JHEP 0612, 003 (2006). W. Altmannshofer, A. J. Buras, S. Gori, P. Paradisi and 
D. M. Straub, Nucl. Phys. B 830, 17 (2010). 

[10] V.M. Abazov et al. (DO Collaboration), Nucl. Instrum. Methods Phys. Res. 
A 565, 463 (2006). 

[11] V.A. Abazov et al. (DO Collaboration), Phys. Rev. D. 82, 032001 (2010). 

[12] V.M. Abazov, et al. (DO Collaboration), Phys. Rev. D 74, 092001 (2006). 

[13] B. Aubert et al. (Babar Collaboratio n), Phys. Rev. Lett. 96, 251802 (2006); 
B. Aubert et al. (Babar Collaboration), arXiv:hep-ex/0607091 (2006). 

[14] E. Nakano et al. (Belle Collaboration), Phys. Rev. D 73, 112002 (2006). 

[15] E. Barberio et al. (HFAG). la7Xiv:0808. 12971 [hep-ex] (2008). 

[16] V.M. Abazov et al. (DO Collaboration), Phys. Rev. D 82, 012003 (2010). 



10 



